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Abstract 

Elliptical rotation is the motion of a point on an ellipse through some angle about a 
vector. The purpose of this paper is to examine the generation of elliptical rotations and to 
interpret the motion of a point on an elipsoid using elliptic inner product and elliptic vector 
product. To generate an elliptical rotation matrix, first we define an elliptical ortogonal 
matrix and an elliptical skew symmetric matrix using the associated inner product. Then we 
use elliptic versions of the famous Rodrigues, Cayley, and Householder methods to construct 
an elliptical rotation matrix. Finally, we define elliptic quaternions and generate an elliptical 
rotation matrix using those quaternions. Each method is proven and is provided with several 
numerical examples. 

Keywords : Elliptic Quaternion, Rotation Matrices, Rodrigues Formula, Cayley Transfor¬ 
mation, Householder Transformation. 

MSC 2010: 15A63, 15A66, 70B05 , 70B10, 70E17, 53A17. 

1 Introduction 

A rotation is an example of an isometry, a map that moves points without changing the 
distances between them. A rotation is a linear transformation that describes the motion of a rigid 
body around a fixed point or an axis and can be expressed with an orthonormal matrix which 
is called a rotation matrix, n x n rotation matrices form a special orthogonal group, denoted 
by SO (n), which, for n > 2, is non-abelian. The group of n x n rotation matrices is isomorphic 
to the group of rotations in an n dimensional space. This means that multiplication of rotation 
matrices corresponds to composition of rotations. Rotation matrices are used extensively for 
computations in geometry, kinematics, physics, computer graphics, animations, and optimization 
problems involving the estimation of rigid body transformations. For this reason, the generation 
of a rotation matrix is considered to be an important problem in mathematics. 

In the two dimensional Euclidean space, a rotation matrix can easily be generated using 
basic linear algebra or complex numbers. Similarly, in the Lorentzian plane, a rotation matrix 
can be generated by double (hyperbolic) numbers. In higher dimensional spaces, obtaining a 
rotation matrix using the inner product is impractical since each column and row of a rotation 
matrix must be a unit vector perpendicular to all other columns and rows, respectively. These 
constraints make it difficult to construct a rotation matrix using the inner product. Instead, 
in higher dimensional spaces, rotation matrices can be generated using various other methods 
such as unit quaternions, the Rodrigues formula, the Cayley formula, and the Householder 
transformation. We will give a brief review of these methods and use elliptical versions of these 
methods later in the paper. 

1. A unit quaternion : Each unit quaternion represents a rotation in the Euclidean 3-space. 
That is, only four numbers are enough to construct a rotation matrix, the only constraint being 
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that the norm of the quaternion is equal to 1. Also, in this method, the rotation angle and 
the rotation axis can be determined easily. However, this method is only valid in the three 
dimensional spaces (0, m- In the Lorentzian space, timelike split quaternions are used 
instead of ordinary usual quaternions m, m)- 

2. Rodrigues Formula : An orthonormal matrix can be obtained using the matrix exponential 
e eA where A is a skew symmetric matrix and 8 is the rotation angle. In this method, only three 
numbers are needed to construct a rotation matrix in the Euclidean 3-space m, m, Esi and, 
[29]).The vector set up with these three numbers gives the rotation axis. This method can be 
extended to the n dimensional Euclidean and Lorentzian spaces (Isa, is], 0 , m and, m)- 

3. Cayley Formula : The formula C = (/ + A) (/ — A) 1 gives a rotation matrix, where A is 
a skew symmetric matrix. Rotation matrices can be given by the Cayley formula without using 
trigonometric functions. The Cayley formula is an easy method but it doesn’t give the rotation 

angle directly (nsi, m, Isa, 0 and, [32]). 

4. Householder Transformation : The Householder transformation gives us a reflection 
matrix. We can obtain a rotation matrix using two Householder transformations. This method 
is an elegant method but it can be long and tedious. Also, the rotation angle has to be calculated 
separately. This transformation can be used in several scalar product spaces m , eq], m, m 
and, [H]). 

Details about generating rotation matrices, particularly in the Euclidean and Lorentzian 
spaces, using these methods can be found in various papers, some of which are given in the 
reference section. Those authors mostly studied the rotation matrices in the positive definite 
scalar product space whose associated matrices are diag(±l, • • • , ±1), and interpreted the results 
geometrically. For example, quaternions and timelike split quaternions were used to generate 
rotation matrices in the three dimensional Euclidean and Lorentzian spaces where the associated 
matrices were diag (1,1,1) and diag (—1,1,1), respectively. In these spaces, rotations occur on 
the sphere x 2 + y 2 + z 2 = r 2 or the hyperboloids — x 2 + y 2 + z 2 = ±r 2 . That is, Euclidean and 
Lorentzian rotation matrices help us to understand spherical and hyperbolic rotations. In the 
Euclidean space, a rotation matrix rotates a point or a rigid body through a circular angle about 
an axis. That is, the motion happens on a circle. Similarly, in the Lorentzian space, a rotation 
matrix rotates a point through an angle about an axis circularly or hyperbolically depending on 
whether the rotation axis is timelike or spacelike, respectively. 

In this paper, we investigate elliptical rotation matrices, which are orthogonal matrices 
in the scalar product space, whose associated matrix is diag(ai, 02 , 03 ) with 01 , 02,03 E M + . 
First, we choose a proper scalar product to the given ellipse (or ellipsoid) such that this ellipse 
(or ellipsoid) is equivalent to a circle (or sphere) for the scalar product space. That is, the 
scalar product doesn’t change the distance between any point on the ellipse (or ellipsoid) and 
origin. Interpreting a motion on an ellipsoid is an important concept since planets usually 
have ellipsoidal shapes and elliptical orbits. The geometry of ellipsoid can be examined using 
affine transformations, because of an ellipsoid can be considered as an affine map of the unit 
sphere. For example, for the ellipsoid £ 2 = |x G K 3 : x*Ax < 1} and the unit sphere S 2 = 
|x G K 3 : ||x|| = x 4 x < l} , we can write £ 2 = T (S' 2 ) using the affine transformation T(x) = 
Ax + c, x E E 2 . Then we have, Vol(£ 2 ) = Vol(T(S 2 )) = -^det QVol (S' 2 ) = v / detQ 47 r /3 where 
Q = A Ah 

The aim of this study is to explain the motion on the ellipsoid 

x 2 y 2 z 2 

-h —-1-= 1 

a 2 b 2 c 2 ’ 

as a rotation, using the proper inner product, vector product and elliptical orthogonal matrices. 
In this method, the elliptical inner product, the vector product and the angles are compatible 
with the parameters 8 and (3 of the parametrization ip(9,/3) = (a cos 8 cos v, b cos 8 sin j3, c sin 8). 
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We use the classical methods to generate elliptical rotation matrices. In the Preliminaries 
section, first we explain how to define a suitable scalar product and a vector product for a 
given ellipsoid. Then we introduce the symmetric, skew symmetric and orthogonal matrices 
in this elliptical scalar product space. Finally, we examine the motion on an ellipsoid using 
elliptical rotation matrices. In section 3, we generate the elliptical rotation matrices using various 
classical methods (such as, Cayley formula, Rodrigues formula and Householder transformation) 
compatible with the defined scalar product. Furthermore, we defined the elliptic quaternions 
and generate elliptical rotations using unit elliptic quaternions. 


2 Preliminaries 


We begin with a brief review of scalar products. More informations can be found in ([!), 
E! and, 0). Consider the map 

(u, v)-X H(u,v) 


for u, v £ M n . If such a map is linear in each argument, that is, 

B (au + fev, w) = aB (u, w) + bB (v, w), 

B (u, cv + dw) = cB (u, v) + dB (u, w), 

where, a, b, c, d £ M and u, v, w £ M n , then it is called a bilinear form. Given a bilinear form on 
M n , there exists a unique H £ M nxn square matrix such that for all u, v £ M n , B (u, v) = iTflv. 
H is called ” the matrix associated with the form” with respect to the standard basis and we will 
denote B (u, v) as Bq (u, v) as needed. A bilinear form is said to be symmetric or skew symmetric 
if B (u, v) =B( u,v) or B (u, v) = — B (u, v), respectively. Hence, the matrix associated with a 
symmetric bilinear form is symmetric, and similarly, the associated matrix of a skew symmetric 
bilinear form is skew symmetric. Also, a bilinear form is nondegenerate if its associated matrix 
is non-singular. That is, for all u £ M n , there exists v £ M n , such that B (u, v) ^ 0. A real scalar 
product is a non-degenerate bilinear form. The space equipped with a fixed scalar product 
is said to be a real scalar product space. Also, some scalar products, like the dot product, have 
positive definitely property. That is, B (u, u) > 0 and B (u, u) = 0 if and only if u = 0. Now, 
we will define a positive definite scalar product, which we call the H-inner product or elliptical 
inner product. 

Let u = (ui,U 2 ,..., u n ), w = (wi,W 2 , •••, w n ) £ M n and a\, 02 ,.., a n £ M + . Then the map 
6 : M n x M" £ M, B (u, w) = a\u\w\ + 02 ^ 2^2 + • • • + a n u n w n 


is a positive definite scalar product. We call it elliptical inner product or H-inner product. The 
real vector space M n equipped with the elliptical inner product will be represented by 02 . an - 
Note that the scalar product B (u, v) can be written as B (u, w) = u f Ow where associated matrix 
is 


n = 


a\ 0 
0 a 2 


0 

0 

0 


( 1 ) 


0 0 ••• a n 


The number \J det fl will be called ” constant of the scalar product ” and denoted by A in the rest 
of the paper. The norm of a vector associated with the scalar product B is defined as ||u|| B = 
V® (u, u). Two vectors u and w are called H-orthogonal or elliptically orthogonal vectors if 
B (u, w) = 0. In addition, if their norms are 1, then they are called H-orthonormal vectors. 
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If {ui, U 2 ,u n } is an ,6-orthonormal base of R” a2 an , then det (ui, U 2 ,u n ) = A . The 
cosine of the angle between two vectors u and w is defined as, 


cos# = 


6 (u, w) 

u ll b II w IIb 


where 0 is compatible with the parameters of the angular parametric equations of ellipse or 
ellipsoid. 


Let B be a non degenerate scalar product, 17 the associated matrix of 6, and R £ M nxn is 
any matrix. 

i) If B (Ru, i?w) = 6 (u, w) for all vectors u, w £ M n , then R is called a 6-orthogonal 
matrix. It means that orthogonal matrices preserve the norm of vectors and satisfy the matrix 
equality R t flR = 17. Also, all rows (or columns) are 6-orthogonal to each other. We denote the 
set of 6-orthogonal matrices by Og (n). That is, 

Og (n) = {R £ M nxn : R^R = Q and det R = ±1}. 

Ok (n) is a subgroup of Gig (n). It is sometimes called the isometry group of M n associated 
with scalar product B. The determinant of a 6-orthogonal matrix can be either —1 or 1. If 
det R = 1, then we call it a 6-rotation matrix or an elliptical rotation matrix. If det R = —1, we 
call it an elliptical reflection matrix. Although the set Og (n) is not a linear subspace of M nxn , 
it is a Lie group. The isometry group for the bilinear or sesquilinear forms can be found in [I]. 
The set of the 6-rotation matrices of can be expressed as follows: 

SOg (n) = {Re M nxn : R^R = 17 and det R = 1}. 


SOg (n) is a subgroup of Og (n). 

ii) If 6 (Shi, w) = 6 (u, Sw) for all vectors u, w £ M n , then S is called a 6-symmetric 
matrix. It satisfies S t i 7 = 17 S. The set of 6-symmetric matrices, defined by 


I = {S £ M nxn : 6 (Shi, w) = 6 (u, Sw) for all u, w £ IT} 

is a Jordan algebra |1]. It is a subspace of the vector space of real nxn matrices, with dimension 
n (n + 1) /2. Any 6-symmetric matrix in a „ can be defined as 


S 


A a, ; 


( 2 ) 


where = aji and a t j £ M. 

iii) If 6 (Tu, w) = — 6 (u, Tw) for all vectors u, w £ M n , then T is called a 6-skew- 
symmetric matrix. Also, T 4 I7 = —17 T. The set of 6-skew symmetric matrices, defined by 


L = {T £ M nxn : 6 (Tu, w) = -6 (u, Tw) for all u, w £ V} 


is a Lie algebra [lj. It is a subspace of the vector space of real nxn matrices, with dimension 
n(n — l)/2, as well. Any 6-skew-symmetric matrix in M” a2 an can be defined as, 


T = Mnxn with tij = 


An 


v 


a 

—An 


ij 


ai 

0 


i > j 

i < j 
i = j 


(3) 


where = aji and aij £ M. 


For example, in the scalar product space Mj* a Q , the symmetric and skew symmetric matrices 
are 
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S = A 


an/ai 

x/a\ 

y/o\ 


0 

x/a\ 

y/ai 

x/a 2 

< 122 / 0-2 

z/a 2 

and T = A 

—x/a 2 

0 

z/a 2 

y/a 3 

z/a 3 

033/03 _ 


_ -y/a 3 

~z/o3 

0 


Note that, even if we omit the scalar product constant A in S or T, they will still be symmetric or 
skew symmetric matrix, respectively. But then, we cannot generate elliptical rotation matrices 
using the Rodrigues and Cayley formulas. So, we will keep the constant A. 

Now, we define the elliptical vector product, which is related to elliptical inner product. 
Let Ui = (tin, Ui 2 , ..., Ui n ) E M n for i = 1,2,..., n — 1 and ei, e2, ..., e n be standard unit vectors 


for B. Then, the elliptical vector product in 


„ is defined as, 

(in 7 


x M n x • • • x M n 


•5 CLn ? 


(ui,u 2 , ...,u n ) 

V (ui X U2 X U3 X ■ 

■ X u n _i) 


ei/ai 

^2/02 

e 3/03 

®n/ 



Uu 

Ul2 

U13 

^1 n 


V (ui x U2 x U3 x • ■ ■ x u n _i) = A det 

u 2 i 

U 2 2 

U23 

'U'2 n 

(4) 


^(n—1)1 

u (n- 1)2 

^(n— 1)3 

^{n— l)n 



The vector V (ui x 112 x U3 x ■ ■ ■ X u n _i) is ^-orthogonal to each of the vectors ui, U2, U3, ..., u re _i. 


3 Generating an Elliptical Rotation Matrix 

In this section, we will generate elliptical rotation matrices using the elliptical versions of 
the classical methods. For a given ellipse in the form 

(E) : \a\x 2 + \ci2y 2 = 1 , A, ai, cl 2 £ R + ( 5 ) 

we will use the scalar product B (u, w) = a\U\W\ + a 2 u 2 W2- That is, our scalar product space is 
R 2 a . An elliptical rotation matrix represents a rotation on (E) or any ellipse similar to (E). 
Recall that ellipses with the same eccentricity are called similar. Since the shape of an ellipse 
depends only on the ratio 01 / 02 , A in Equation ([5]), does not affect the rotation matrix. 

3.1 Rodrigues Rotation Formula 

The Rodrigues rotation formula is a useful method for generating rotation matrices. Given a 
rotation axis and an angle, we can readily generate a rotation matrix using this method. SO(n) 
is a Lie group and the space of skew-symmetric matrices of dimension n is the Lie algebra 
of SO(n). We denote this Lie algebra by so(n). The exponential map defined by the standard 
matrix exponential series e A connects it to the Lie group. For any skew-symmetric matrix A, the 
matrix exponential e A always gives a rotation matrix. This method is known as the Rodrigues 
formula. 


3 . 1.1 Elliptical Rotations In the Plane 


According to 0 a skew symmetric matrix can be expressed as 


T = 


0 —yfwil ^Jai 

y/al/y/a^ 0 


The equality = —QT is satisfied. The characteristic polynomial of T is, P (x) = x 2 + 1. So, 
T 2 = —I. We can obtain the elliptical rotation matrix easily using the matrix exponential. 
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Theorem 1 Let T be a B-skew symmetric matrix. Then, the matrix exponential 


R$ = e Te = I + (sin 6) T + (1 - cos 9) T 2 = 


cos 9 


y/^2 . Q 

-sin 0 

y/ai 


LL. sin 6 cos 9 

L \/®2 


gives an elliptical rotation along the ellipse \a\x 2 + Xa 2 y 2 = 1, A, a\,a 2 G R + . That is, Rq is 
an elliptical rotation matrix in the space a . 

Remark 1 All similar ellipses have identical elliptical rotation matrices. 

X 2 

Example 1 Let's consider the ellipse (Ei) : — + — = 1 with the parametrization a (6) = 

9 4 

(3 cos 9, 2 sin 9). Let’s take the points 


A = a (vr/4) = 


,V2] and B = a (tt/4 + tt/3) = (3\/2/4 - 3^6/4, ^2/2 + V6/2). 


That is, if we rotate the point A elliptically through the angle n/S, we get B. Now, let’s find 
elliptical rotation matrix for these ellipses using the above theorem and get same results. To 
calculate the elliptical rotation matrix, first, we choose the elliptical inner product 


B (u, w) 


U\W\ U2VJ2 

9 + 4 


in accordance with the ellipse (Ei) such that u = (u\,U 2 ) and w = (w\,W 2 ). That is, our space 
is M^ g 1 , 4 , A = 1/6 and the B-skew symmetric matrix has the form 


0 -3/2 

2/3 0 


Note that, T 2 = —L. We can obtain elliptical rotation matrix as, 


r>B _ JT 
Rq - e 


cos 9 
2 sin 9 
3 


3sin0 

2 

cos 9 


where Rq is a B-orthogonal matrix in ^1/9,1/4- Namely, the equalities det R® = 1 and (RqY (i?f) 
fl are satisfied. For 9 = 7r/3, we get 


R 


B 

7r/3 


1/2 —3\/3/4 

V3/3 1/2 


So, if we rotate the point A elliptically, we get B = Rq (A) = (3\/2/4 — 3\/6/4, y/2/2 + \/6/2). 
Thus, we get same result using elliptical rotation matrix for (Ei ). Note that, \\Rq (A)|| b = 1 
and the angle between x =OA and y =oY is 


cos U = 


B(x,y) 


b llylls 


X 


R can be seen that the elliptical rotation matrix Rq can be also used to interpret the motion on 


x 2 y 2 

a similar ellipse -1-= 1 to (Ei). 

36 16 
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Figure 1: a (#) = (3 cos #, 2 sin #) and (3(9) = (6 cos #,4 sin #) 


3.1.2 3-Dimensional Elliptical Rotations 

Let’s take the ellipsoid a\x 2 + a 2 y 2 + a^z 2 = 1. The scalar product for this ellipsoid is 


B (u, w) = aiUlWl + CL2U2W2 + CI3U1W3, 
for u = (m,U2, U3 ) and v = (ui, V2, V3) . Also, the vector product is 



ei/ai e 2 /a2 e 3 /a3 


0 — u 3 /a\ 112/ai 


" Vl 

V (u x v)=A det 

Ui U 2 U3 

=A 

143/02 0 -ixi/a 2 


V2 


V\ v 2 V 3 


_ -ix 2 /a 3 ui/a 3 0 


. v 3 _ 


where A = sJTFfafaf. The matrix 


T = A 


0 -u 3 /ai 112/(11 
U3/02 0 ~ ui / a 2 

-112/0,3 u\/ o 3 0 


( 6 ) 


is skew symmetric in . That is, T f i\ = —LIT. So, the vector product in a3 can 

be viewed as a linear transformation, which corresponds to multiplication by a skew symmetric 
matrix. The characteristic polynomial of T is, P (x) = x 3 + || u|| 2 x whose eigenvalues are x\ = 0 
and X 2,3 = ± ||u|| i. According to characteristic polynomial T 3 + 11u11 2 T = 0. So, if we take a 
unit vector u £ K 3 Q Q , we get T 3 = — T and we can use Rodrigues and Cayley formulas. 


Theorem 2 Let T be a skew symmetric matrix in the form © such that u = ( 1 x 1 , 1 x 2 , 143 ) £ 
M 3 i a2 a3 is a unit vector. Then, the matrix exponential 

R®’ u = e Te = I + (sin #) T + (1 - cos 6) T 2 

gives an elliptical rotation on the ellipsoid a\x 2 + a^y 2 + a 3 Z 2 = 1. Furthermore, the matrix Rq 
can be expressed as 

a\u\ + (l — a\uf) cos 9 — Au ^ in6 — 021x11x2 (cos# — 1) Au ^ m9 — 031x11x3 (cos# — 1) 

Au ^ sm6 — a\U\U2 (cos# — 1) 02^2 + (l — a 2 U 2 ) cos# — AM l^ ln6> — 031X21X3 (cos# — 1) 

— Sin6> _ aiUlU3 (cos# — 1 ) A “i 3 Bln " g — 02 U 2 U 3 (cOS # — 1 ) 031X3 + (l — 031X3) COS# 

( 7 ) 

where u = (ixi, 112 , 113 ) is the rotation axis and # is the elliptical rotation angle. 
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Proof. Since u is a unit vector in 


~ai,a2,<X3 


, we have T 3 = —T. So, we get 


Bu or „ 9 2 T 2 -9 3 T -9 4 T 2 9 3 T 9 6 T 2 

Re' = e 9T = I+ 9T+ + —— + — + 


2 ! 3 ! 

( 9 3 9 5 

= / + r ( 9 -W + ^-' 

= I + t(9- 9 ^ + — -^j+T 2 

v 3! 5! ) 

= 1 + (sin0)T + (1 -cos0)T 2 . 


4! 5! 

, 9 2 0 4 0 6 

+T~ - 1 - 

1 2! 4! 6! 


6 ! 


, . 1 9 2 9 4 0 6 

1 ” 1 1 " 2! + 4! " 61 + 


If we expand this formula using —C 12 U 2 — a%u 2 = aiu 2 — 1, we can obtain the rotation matrix as 

0 . - 



Figure 2: 


x 2 y 2 z 2 

Example 2 A parametrization of the ellipsoid —H- —j— H- = 1 is 

a (9 , /3) = (2 cos 9 cos /?, 2 cos 9 sin /3,3 sin 9) 


where 9 G [0, 7 r) and (3 G [0, 27t). Let’s take the points 

A=a (30°, 30°) = ^3/2, \/3/ 2 , 3/2^ and B=a (120°, 30°) = (-\/3/2,- 1 / 2 , 3\/3/2) 

on t/ie ellipsoid. Let’s find the rotation matrix which is rotate the point A to B elliptically. We 
have ai = a 2 = 1/4 and 03 = 1/9. £ 0 , A = 1/12. First, using the vector product of x = OA and 
y = o£ in M 3 /4 1 / a 4 ^, we find the rotation axis u. 


4i 4j 9/e 

1 3/2 \/ 3/2 3/2 

' J -\/ 3/2 3 i/ 3/2 


(l, -x/3, 0 ) . 


.Since V (x x y) is unit vector in M 3 y 4 4 ^ g , we get u = ( 1 , — \/3,0). Thus, we obtain the elliptical 
rotation matrix 


9 cos 9 + 3 
3\/3 (cos 9 — 1) 
9y / 3sin0 


3\/3 (cos 9 — 1) 
3 cos 0 + 9 
9 sin 9 


—4^3 sin 9 
—4 sin 9 
12 cos 9 


ReT(9) 


1 

12 























by using &■ This matrix describes an elliptical rotation on a great ellipse such that it is 
intersection of the ellipsoid and the plane passing through the origin and B-orthogonal to u. It 
can be easily found that equation of the plane is x = V^y. So, R “ represents an elliptical rotation 


1 


over the the great ellipse is y 2 + -z 2 = 1, y = V?>x. 
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Also, the elliptical rotation angle is 7r/2, 

since cos 6 = B (x, y) = 0 (Figure 2b). Thus, we find 

3 —3\/3 


R u _ 1 

K/2 - 


—3\/3 
9\/3 


9 

9 


-4\/3 

-4 

0 


( 8 ) 


The matrix (8), rotates the point A to the point B elliptically over the great ellipse — + — = 1, 
_ 2 9 

y = x. 

Remark 2 The eigenvalues of the matrix ^ are, x\ = e l6 , X 2 = e~ ld and x 3 = 1. Also, the 
eigenvector corresponding to 1 is the vector u, the rotation axis of the motion. 


3.2 Elliptical Cayley Rotation Matrix 

In 1846, Arthur Cayley discovered a formula to express special orthogonal matrices using 
skew-symmetric matrices. It is called Cayley rotation matrix. In this section, we will describe 
the Cayley rotation matrix for any ellipsoid whose scalar product is B. We call it £>-Cayley 
rotation matrix. 

The Cayley rotation matrix is a useful tool that gives a parameterization for rotation matri¬ 
ces without the need to use trigonometric functions. Let A be an n x n skew symmetric matrix 
where (/ — A) is invertible. Cayley formula transforms the matrix A into (I + A) (I — A)~ l . 
SO (n) is isomorphic to so(n) via the Cayley formula where so (n) denotes the space of skew 
symmetric matrices, usually associated with the Lie algebra of the transformation group defined 
by rotations in SO(n). That is, the Cayley map is defined as 

C : so (n) —> SO (n) 

A -» C (A) = (I + A) (I - A)- 1 . 


Note that, if a matrix M has the eigenvalue —1, then (/ — M) is not invertible. But, since 
A is skew symmetric, all its eigenvalues are purely imaginary in the Euclidean space and (/ — A) 
is invertible. That is, the Cayley formula is well-defined for all skew symmetric matrices. The 
inverse of the Cayley map is given by 

C~ l = {I — A) (I + A)- 1 . 


In the Minkowski space, an eigenvalue of a skew symmetric matrix can be —1. In this case, 
the Cayley formula is not valid. For detailed information on the Cayley formula in Minkowski 
space see [18]. 


Theorem 3 Let T be a skew symmetric matrix in the form & such that u = (ui,U 2 ,U 3 ) £ 

^,02,03 • Then > 

R s ’ u = (/ + T) (I - T)" 1 

is an elliptical rotation matrix on the ellipsoid a\x 2 + 02 y 2 + a^z 2 = 1 where u is the rotation 
axis. Furthermore, the matrix R® ,u can be written in the form 


l 

l+l! u !le 


a\u 2 —u 2 a 2 ~ ^^ 3+1 
^ 3 -+a 1 uiu 2 +uiu 2 a 1 
aruiu 3 -'^^+uiu 3 ai 


a 2 u 1 u 2 - +U 1 U 2 Q 2 

0 0 o 

a2U2—ufa\—u^a 3 +l 

‘^ L +a2U2U 3 +u 2 u 3 a2 


^ ( z +a 3 uiu 3 +uiu 3 a 3 
a 3 U2U 3 —^^ x +U2U 3 a 3 
a 3 u 2 —u 2 ai~ tt 2 a 2 +l 


(9) 
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Proof. Since, T is a R-skew symmetric matrix, we have = —LIT. Also, we can write 
(/ + T) 1 n = n (/ - T) and (/ - T) f Q = Ll(I + T). 


Using these equalities, it can be seen that 

(R B ’ u y n ( R B ’ u ) = ((/ + T) (/ - T)" 1 )* U (I + T) (/ - T)" 1 = n. 

Also, since det (/ + T) = 1 + ||u|| B and det (I — T) _1 = --——, we have detR Bu = 1. That 

i + II u IIb 

is, R b ’ u is an elliptical rotation matrix. The matrix (|9j) can be obtained after some tedious 
computations. ■ 


Remark 


3 The eigenvalues of the matrix 0 are 


Ai — 


1 - 11 u 11 
1 + Hull 


+ 


2 ||u| 


1 + ||u 


12 ’ 


^2 — 


1 - ||u| 
1 + ||u| 


2 ||u|| B i 

1 . 1 p" 

1+ l|u|| 6 


and A 3 = 1 . 


Also, the eigenvector corresponding to 1 is u, which is the rotation axis. The matrix 0 rotates 
a vector elliptically on the ellipsoid a\x 2 -Ua^y 2 + asz 2 = 1 about the axis u = (tti, U 2 , U 3 ) through 
the elliptical angle 9 where 


tan 9 


2 ll u llg 

1 11 112 1 

1 — ll u lls 


( 10 ) 


Example 3 Let’s find the elliptical rotation matrix representing a elliptical rotation on the 
x 2 y 2 

ellipsoid - 1 -f z 2 = 1 about the axis u = (2, 3,1). Using the matrix (9), we find 

4 9 


002 
3/2 0 0 

0 1/3 0 


The elliptical rotation angle corresponding to this matrix, can be found to be —ir/3 from the 
formula (10). 


3.3 Elliptical Householder Transformation 

The Householder transformation was introduced in 1958 by Alston Scott Householder. A 
Householder transformation is a linear transformation in the form 

2 vv 4 

H v (x)=x - 7 —x 

V t V 

where v is a nonzero vector. This transformation describes a reflection about a plane or hyper¬ 
plane passing through the origin and orthogonal to v. Householder transformations on spaces 
with a non-degenerate bilinear or sesquilinear forms are studied in [lj. Every orthogonal trans¬ 
formation is the combination of reflections with respect to hyperplanes. This is known as the 
Cartan-Dieudonne theorem. A constructive proof of the Cartan-Dieudonne theorem for the 
case of generalized scalar product spaces is given by Uhlig HU and Fuller m ■ An alternative 
proof of the Cartan-Dieudonne theorem for generalized real scalar product spaces is given by 
Rodriguez-Andrade and etc. [2\. They used the Clifford algebras to compute the factorization 
of a given orthogonal transformation as a product of reflections. 

Householder transformations are widely used in tridiagonalization of symmetric matrices 
and to perform QR decompositions in numerical linear algebra |3j. Householder Transformation 
is cited in the top 10 algorithms of the 20th century [35]. Generalized Householder matrices 
are the simplest generalized orthogonal matrices. A generalized Householder, or B- Householder, 
matrix has the form 
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H v (x) = x 


2Wn 


for a non-B-isotropic vector v (i.e., v*fiv 0) [13j, [2j. 

In this paper, we use the elliptical version of the Householder transformation to generate an 
elliptical rotation matrix in a real elliptical scalar product space M" a an . Let H v = [hij] nxn 
be a Householder matrix. Then 

, r 2 ViVjQjj 

hij = 6ij - v'Ov 


where 5ij is the Kronecker delta. The Householder matrix H v is a H-synunetric and B-orthogonal 
matrix. To prove these, first we will show that = LlH v : v 4 Hv = £>(v,v) = ||v||g is a 

positive real number. So, 


Kn 


f 2Wn 

y v*Hv 


n = 



2S7vv < H\ 
v*flv ) 


= n 



2vv t H\ 

v*Hv J 


nn v , 


and 7~i v is a B-symmetric. Next, we will show that H 1 v O'Hv = H. 


ninn v 



2 (H v is B-symmetric) 
2vv < H\ / 2Wn\ 

v*Hv ) y v^v J 


f 2 vv 4 H 2vv*H 4 (vv 4 !!) (vv*f2) \ 
y v*Hv v 4 flv + ( v ( Ov) 2 J 

= a 


So, H v is B-orthogonal. 


Using the above transformations, we can write the elliptical reflection matrix as 


1 


'Hv = -—Tw —2aiU2Ui 

-2aiu 3 ui 

for the ellipsoid aix 2 + a 2 y 2 + a 3 z 2 = 1. 


— 2a\v\ — 2 a 2 V\V 2 


-2a 3 uiu 3 


|v|| B — 2 a 2 t >2 — 2 a 3 U 2 f 3 


— 2a 2 v 3 V2 


— 2a 3 U3 


(ID 


Corollary 1 The Householder transformation doesn’t change the length of a vector in Ka lja2ja3 - 
It describes an elliptical reflection about a plane passing through the origin B-orthogonal to v. 
Using the matrix 0, we can express elliptical reflections on any ellipsoid. 


Example 4 Let’s find the elliptical reflection matrix which reflects a point elliptically on the 
ellipsoid 2x 2 + 2 y 2 + z 2 = 1, about the plane which is passes through the origin and is B- 
orthogonal to v = (1, 2, 3). Our elliptical inner product is B (x, y ) = 2xi2/2 + 2 x 22/2 + ^32/3- So, 
= 2 • l 2 + 2 • 2 2 + 1 • 3 2 = 19 and we get 


i2 

11 B 


H 


V 


15/19 -8/19 -6/19 
-8/19 3/19 -12/19 
-12/19 -24/19 1/19 


The determinant of this matrix is — 1 and the equality ’H t v O'H v = H holds. We find the reflection 
of A (1/2,1/2,0) to be B(7/38,-5/38,-18/19). Let’s check it : The middle point of [AB\ is 
(7(13/38,7/38,-9/19). The equation of the plane is B(v,x) = 0 where x = ( x,y,z ). In our 
case, it is 2x + Ay + 3z = 0, so C lies on this plane. Also, it can be easily seen that .4 fl is 
B-orthogonal to the plane. 
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Remark 4 For the n dimensional ellipsoid associated with the scalar product B, the elliptical 
reflection matrix FL V is in the form 


qjB _ 1 

II 112 

I v b 

||v g — 2a\v\ 
—2a\V2V\ 
—2a\v^v\ 

-2a 2 viV2 

. | . 1 2 o 

V g — 202^2 

-2a2V 3 V2 

-2a 3 viv 3 
—2a 3 V2V 3 
v B — 2a 3^3 '• 

2 d n V\V n 
‘2‘ClriV2'Vn 
2cL n V^V n 

_ —2a\v n vi 

-2 a 2 v n V2 

-2 v n a 3 v 3 

• v||| - 2 a n v\ 

where v f = (ui, t> 2 ,. 

v n ). 





Theorem 5 Let x, y G Mg be any two nonzero vectors where ||cc||g = ||y||g. If v = x — y, then 

„ ,n , \ 2vv f fl 

Hv (x) = * " = y - 

Proof. By direct computation, we find 

v*Hv = B (x — y, x — y) = 2 ||x||g — 2y t Llx 


and 


2vv t Tlx = 2 (x — y) (x — yy Qx = (x — y) ( 2 ||cc||g — 2y t Llx 


Using these equalities, we obtain 


„ , 2WQ (x-y)[2\\x\\ B -2y t nx 

( x ) = * - x = x- 


2 ||cc|| e — 2y t Llx 


= x- x + y = y. 


Theorem 6 Let x, y G Mg 6e any two nonzero vectors where ||cc||g = ||y||g. //v = x + y, then 

U® (x) = -y. 


Proof. By direct computation, we see that 

v f Ov = B (x + y, x + y) = 2 ||®||g + 2y t Llx 


and 

2vv t Llx = 2 (x + y) (x + y) 4 fix = (® + y) ^2 ||cc||g + 2y t Llx s j . 

Thus, we get 

"Hv (*) = x- (® + y) = -y • 

Moreover, we have (—y) = y. 


2yy t Q 


2 y (y*fiy) 


(-2/) = -y + ^771772/ = “2/ + -777-“ = 2/- 


y*Hy 


ytfly 


which means that HyH^-\-y (*) = 2/- Hence, we obtain the elliptical rotation matrix 


RS = RyR^+y 

which rotates x to y elliptically. The rotation axis is the vector V (x x y) and the elliptic rotation 
angle can be found using elliptical inner product. ■ 
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Example 6 Let’s find the elliptic rotation matrix that rotates x = (0,0,5) to y = (2,2,3) on 
the ellipsoid 2x 2 + 2 y 2 + z 2 = 1. Using the matrix (11), we get 


H 


x+y 


4 

-1 

-2 ' 

and R y = ^= 

9 

-16 

-12 ' 

-1 

4 

-2 

-16 

9 

-12 

. - 4 

is 

-4 

-3 

25 

-24 

-24 

7 


R = 'Hy'H X + y = 


1 


4-12 
-14 2 

-4 -4 3 


which satisfies R [x) = y, det R = 1, and R t LlR = 12 where 12 = diag (2, 2,1). ./Vole that, we can 
obtain same matrix using 0- 


3.4 Elliptic Quaternions 


Quaternions were discovered by Sir William R. Hamilton in 1843 and the theory of quater¬ 
nions was expanded to include applications such as rotations in the early 20th century. The 
most important property of the quaternions is that every unit quaternion represents a rotation 
and this plays an important role in the study of rotations in 3-dimensional vector spaces. Using 
unit quaternions is a useful, natural, and elegant way to perceive rotations. Quaternions are 
used especially in computer vision, computer graphics, animation, and kinematics. 

Quaternion algebra HI is an associative, non-commutative division ring with four basic 
elements { 1, i, j, k} satisfying the equalities i 2 = j 2 = k 2 = ijk = —1. We can express any 
quaternion q as q = (qi,q 2 ,Q 3 , ( u) = <Zi +< 72 i + <Z 3 j + g 4 k or q = S q + Vg where the symbols S q = q\ 
and Vg = q 2 \ + q?fii Tgqk denote the scalar and vector parts of q, respectively. If S q = 0 then q is 
called a pure quaternion. The conjugate of q is denoted by q, and defined as q = S q — ~V q . The 
norm of a quaternion q = qo + qfi + g 2 j + (/3 k is defined by fiqq = fiqq = \Jq^ + qf + + q 2 and 

is denoted by N q . We say that qo = q/N q is unit quaternion if q 0. The set of unit quaternions 
are denoted by Hi. Every unit quaternion can be written in the form qo = cos 9 + e 0 sin 9 where 
£0 is a unit vector satisfying the equality £q = —l.It is called the axis of the quaternion 0, El- 

Let q and r be two quaternions. Then, the linear transformation R q : H —> H defined by 
Rq ( r ) = q r q l is a quaternion that has the same norm and scalar as r. Since the scalar part of 
the quaternion r doesn’t change under R q , we will only examine how its vector part V r changes 
under the transformation R q . We can interpret the rotation of a vector in the Euclidean 3-space 
using the quaternion product q\ r q~ 1 . 

If q = qo + qfii + (72 j + Q 3 k = cos 9 + £osin0 is a unit quaternion, then, using the linear 
transformation R q (V r ) = q^rQ. 1 , the corresponding rotation matrix can be found as 


Rq 


<?o + Qi - 02 - <?3 -2<7o<?3 + 2gi g 2 ‘i-WU + 2gig 3 

2gig 2 + 2g 3 g 0 Qo ~ Qi + Q 2 ~ Q 3 2 g 2 g 3 - 2gig 0 
2qiq 3 - 2q 2 q 0 Zqiqo + 2g 2 g 3 ql - q 2 - q 2 + q 2 


(4) 


This rotation matrix represents a rotation through angle 29 about the axis £ = (gi,g 2 ,g 3 ) • In 
the Lorentzian space, the rotation matrix corresponding to a timelike quaternion q = qo + q\\ + 
qfii + ©k is, 


Rq 


qi + q z 2 + qi + < 
2g 2 <23 + 2 <24<2 i 
2g 2 <24 - 2 <23<2 i 


2q\q4 ~ 2 q 2 q 3 

q‘i - q‘i ~ qi + ql 

2q 2 qi - 2 g 3 g 4 


-2qiq 3 - 2 q 2 q 4 
-2g 3 <24 - 2<2 2 <2 i 

ql -ql + qi-ti 


where the set of unit timelike quaternions satisfies the properties qq = 1, i 2 = — 1, j 2 = k 2 = 
ijk = 1 (see [16]). 
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3.4.1 Elliptic Quaternions 


To get an elliptical rotation matrix, firs we define the set of elliptic quaternions suitable 
for the ellipsoid a\x 2 + a 2 y 2 + a 3 z 2 = 1. Let’s take four basic elements {1, i, j, k} satisfying the 
equalities 


• 2 *2 i 2 

i = -a 1 , y = — a 2 , k- = — a 3 


and 


ij = 



jk 



kj, 


where 01 , 02,03 G M + and A = s /a\a 2 a 3 . 


ki = 



-ik 


The set of elliptic quaternions will be denoted by H ai)a2!a3 . This set is an associative, non- 
commutative division ring with our basic elements { 1 , i, j , k} . If we take 01 = a 2 = 03 = 1, we 
get the usual quaternion algebra. The elliptic quaternion product table is given below. 



1 

i 

j 

k 

1 

1 

i 

j 

k 

i 

i 

Ol 

Ak/a 3 

-Aj/a 2 

j 

j 

—Ak/a 3 

—a 2 

Ai/ai 

k 

k 

Aj/a 2 

—Ai/ai 

-03 


The elliptic quaternion product of two quaternions p = po + Pii + p 2 j + p 3 k and q = 
1 o + qii + q 2 j + gqk is defined as 

POlO — B (Vp, V 9 ) + PqV q + qoV p + V (Vp X V 9 ) (12) 

where B(V p ,\ q ) and V (V p x V q ) are the elliptical scalar product and the elliptical vector 
product, respectively. If p and q are pure, then 

VI = (Vp, Vq) + V (Vp x V 9 ) 


= - (ciipiqi + ( 12 P 2 Q 2 + a 3 p 3 q3 ) + A 
The elliptic quaternion product for H ai)a2ia3 can be expressed as 


i /01 

j '/«2 

k/a 3 

Pi 

P2 

P3 


12 

13 


pq = 


Po 

-aipi 

—C12P2 

P3A 

- a 3 p 3 

P 2 A 

Pi 

0 <1 

ai 

O' 1 

Pi A 

P2 

<1 

CM CM 

0 <1 
^ £ 

03 

P3 

0.3 

03 

Po 


% 

Qi 

Q2 

Q3 


For example, let p, q G H 2 , 2 ,i- Then, the elliptic quaternion product of p and q defined is 


pq = 


po —2pi -2 p 2 ~P3 

Pi Po ~P3 P2 

P2 P3 Po -Pi 

P3 -2 P2 2pi po 



’ 10 ’ 


11 


12 


. 13 . 
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For p = 1 + 2i + 3j + 4k and q = 2 + 4i + j + 3k, we get pq = (—32,13,17, —9). This can also 
be calculated using the product table 



1 

i 

j 

k 

1 

1 

i 

j 

k 

i 

i 

-2 

2 k 

-j 

j 

j 

-2k 

-2 

i 

k 

k 

j 

— i 

-1 


Remember that the algebra is formed by a vector space V equipped with a quadratic form Q 
with the following equalities 


v 2 = Q (v); 

uv + vu = 2 Bq (u, v) 

is called a Clifford algebra and is denoted by Cl(Y,Q). If {ei,...,e n } is a base for an n- 
dimensional vector space V, then Cl (V, Q ) is formed by the multivectors 

{l}U{e il e ia ...e ifc :l < ii < ... < < n, 1 < k < n} 

with dim (Cl (V, Q)) = 2 n . Since the Clifford product of two even multivectors is an even mul¬ 
tivector, they define an even subalgebra of Cl (V, Q). The even subalgebra of an n-dimensional 
Clifford algebra is isomorphic to a Clifford algebra of (n — 1) dimensions and it is denoted 
by Cl + (V, Q) . The Hamiltonian quaternion algebra HI is isomorphic with the even subalgebra 
Clto = Cl (M 3 , Q = xI + x 2 + £ 3 ) by {1, e 2 e 3 — > j, eie 3 —> k, eie 2 —> i} and the split quater¬ 
nion algebra HI is isomorphic with the even subalgebra Cl\ x = Cl (M 3 ,Q = — x\ + x\ + # 3 ) 
by {I,e 2 e 3 —> i,e 3 ei —> k, eie 2 —> j} 0. Similarly, the elliptic quaternion algebra is an even 
subalgebra of the Clifford algebra 

Cl (M 3 ) = {q = qo + eiqi + e 2 q 2 + e 3 q 3 : e\ = ai, e| = a 2 , e| = a 3 , e t e 3 + = 0 } 

associated with the nondegenerate quadratic form Q (x) = a\x\ + + a 3 x\ and is de¬ 

noted by Cl + (M 3 ,aixf + 02^2 + «3®|) , or shortly C£ + (M 3 , a2 a ,J . H aija2ja3 is isomorphic to 

Cl + (K 3 lja2ja3 ) with { 1 , ^e 2 e 3 -> i, ^eie 3 -> j, ^eie 2 -> k}. 

For the quadratic form Q = a\x\ + aix| + a\x 3 , recall that the elliptical inner product can be 
obtained by using the equality 

>3q (x, y) = ^ [Q (x + y) - Q (x) - Q (y)]. 

So, we get Bq (x, y) = aiaqyi + a 2 x 2 y2 + a 3 x 3 y 3 for x = (xi,x 2 , x 3 ) and y = (yi, y 2 , y 3 ) ■ Thus, 
we can construct an elliptic quaternion algebra for any elliptical inner product space. 
Conjugate, norm and inverse of an elliptic quaternion q = qo + qfi + q 2 j + q 3 \i can be defined 
similar to usual quaternions : 

q = (jo- qii- N q = ^/qq = sffq = \Jq% + aiqf + a 2 ql + a 3 q f, q~ l = J-. 

Also, each elliptic quaternion q = qo + qii + q 2 j + ( 73 k can be written in the form 

q 0 = N q (cos 0 + e 0 sin 0 ) 


where 
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a 90 A ■ a V a i9i + a 2Q 2 + “393 

cos u = —— and sm 6 = — - 


N q N q 


Here, £q = 


( 91 , 92 , 93 ) 


is a unit vector in the scalar product space 


ai ,02,03 


satisfying 


■\/“i9? + a 2 gf + a 3 g| 

the equality £g = —1. It is called the axis of the rotation. For example, if q = 1 + 2i + j + 5k G 
H 2 , 2 ,i) then N q = \/l 2 + 2 • 2 2 + 2 • l 2 + 1 • 5 2 = 6 and we can write 

1 , v/35(2,l,5) „ (2,1,5) . 

6 6 ^35 V35 


where £q = 


v/35 


(2,1, 5) is a unit vector in M 3 ,2 i with £q = —1. 


Theorem 7 Fac/i a unit elliptic quaternion represents an elliptical rotation on the an ellipsoid. 

If 

q = qo + 9ii + 92 j + <? 3 k = cos 9 + So sin 9 G H ai)a2)a3 

is a unit elliptic quaternion, then the linear map Rg (v) = qvq~ 1 gives an elliptical rotation 
through the elliptical angle 26, about the axis £q, where v G M 3 . The elliptical rotation matrix to 
corresponding to the quaternion q is 


K = 


9o+“i9i - a 2 ql - a 3 ql 2a 2 qiq 2 ~ 2 

A 2 2 


2a\q\q 2 + 2- 
2a 1 q 1 q 3 - 2 


9o93A 

a i 

.2 „ ,,2 


2a 3 qiq 3 + 2 


02 A 

©92 A 
03 


90 “ 0191+0292 - °393 2 “39293 - 2 


2a 2 9293 + 2 


go 9i A 
03 


9o 92 A 

o 1 

9091A 
02 


9o - “i9i - 0292+0393 


(13) 


Proof. It can be seen that Rg is a linear transformation and preserves the norm. Using the 
equalities, 


Rg (i) = (aiqlql~qla 2 -qla 3 ) i + 2 (aiqiq 2 +qoq 3 J+ 2 ^9i93ai-go92 ^ 

Re (j) = 2 ^029192-9093^°^^ i+ ( 029 ! + 9o-9i“i-g|“3) j+2 ^a 2 9293+9o9i \j 
Rg (k) = 2 ( a 3 q 1 q 3 +q 0 q 2 


a 3 a 2 


k, 

©02 \ j 

03 J 


i+2 (“39293-9091 \j a ~^j j+ (“393 + 9o-9?ai-92a 2 ) k, 


we can obtain (13). So, the rotation matrix (13) is an elliptical rotation matrix on the ellipsoid 
a\x 2 + a 2 y 2 + a 3 z 2 = 1. That is, the equalities det-R# = 1 and RgQRg = H are satisfied. Also, 
note that, if we take a\ = a 2 = a 3 = 1, the standard rotation matrix is obtained. Now, let’s 
choose an orthonormal set {so, £ 1 , e 2 } satisfying the equalities 

V(e 0 xei) = e 2 , V(e 2 xe 0 )=ei, V(£iX£ 2 ) = £o- 

If e is a vector in the plane of the £0 and e \, we can write it as 

£ = cos aeo + sin a£i. 

To compute Rg (e) = 9 £g _1 , let’s find how £0 and £1 change under the transformation Rg. Since 
V 9 is parallel to e 0 , we have qeq = £09 by ([ 4 ]) and R q (e 0 ) = g£o 9 1 = £ 0 99 1 = £o- So, £0 is 
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not changed under the transformation R g . It means that £o is the rotation axis. On the other 
hand, 

R q (si) = qeiq- 1 

= (cos 9 + £q sin 9) e\ (cos 9 — e o sin 9) 

= £i cos 2 6 — cos 9 sin 9 (si£o) + cos 9 sin 9 (£o£i) — (£o£i) £o sin 2 9. 

Since we know that £i£o = V (si x £o) = —V (^o X £i) = — £o£i = —£2 for orthogonal, pure 
quaternions, we obtain 

R q (£i) = £i cos 2 9 + (£i£o) £o sin 2 9 + 2£ 2 cos 9 sin 6 
= £i cos 2 9 + e i£q sin 2 9 + 2e2 cos 9 sin 9 
= £i cos 29 + £2 sin 29 


That is, £ is rotated through the elliptical angle 29 about £q by the transformation R q (£). m 


Corollary 1 All elliptical rotations on an ellipsoid can be represented by elliptic quaternions 
which is defined for that ellipsoid. 


Example 6 Let's find the general elliptical rotation matrix for the ellipsoid 2x 2 + 2y 2 + z 2 = 1. 
Using (13), we obtain, 



Qo + 2 <?i - 2 <?2 “ $ - 2go93 2g 0 <?2 + Zqm 

2<?o<?3 + 4<?i<72 ql - 2q( + 2 q% - q( 2q 2 q 3 - 2q 0 q 1 
4qm - 4<?o<?2 4g 0 gi + 4®® ql - 2 q\ - 2q\ + ql 


Here, det Re = ((/q + 2q\ + 2q 2 + ql) 3 = 1 and R^UlR = 0 where Ll = diag (2, 2,1). For example, 
the unit quaternion q = (0,1/2,1/2, 0) represents an elliptical rotation on the ellipsoid 2x 2 + 
2 y 2 + z 2 = 1 through the elliptical angle it, about the axis (1/2,1/2, 0). And the elliptical rotation 
matrix is 


R 


9 

7T 


0 1 0 

1 0 0 

0 0-1 


3.5 An Algorithm 

Generating 3-dimensional rotation matrix that rotates x = (xi,yi,yi) to y = ( x 2 ,y 2 ,z 2 ) ellipti- 
cally on the ellipsoid a±x 2 + a 2 y 2 + a 3 z 2 = 1 

Step 1. Write, a = (ai, a 2 , 03 ) for the given ellipsoid Qqx 2 + a 2 y 2 + a 3 z 2 = 1 where ai, a 2 , 03 G 
R+. 

Step 2. Define Scalar Product B, norm of a vector and Scalar product constant A as follows : 
B (x,y, a) = B (xi, yi, z 3 , x 2 , y 2 , z 2 , ai, a 2 , a 3 ) = aixi;c 2 + a 2 y\y 2 + a 3 z 3 z 2 , 

A f (x,a) = y/B ( x , x, a), 

A = y/a\a 2 a 3 

where x = (xi,yi,yi) and y = (x 2 , y 2 , z 2 ). 

Step 3. Define Vector Product V as 

V(x,y, a) = V {xi,yi,zi,X 2 ,y 2 ,z 2 ,ai,a 2 ,a 3 ) 

— (yiZ 2 - V 2 Zi ), — (-X 1 Z 2 + x 2 zi ), — (xiy 2 - x 2 yi) 
ai a 2 a 3 
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Step 4. Choose the vectors x = (a;i,yi,yi) and y = ( x 2 ,y 2 ,z 2 ) to find the elliptical rotation 
matrix that rotates x to y elliptically on the ellipsoid. 

Step 5. Find, V (x, y, a) and norm of the vectors x, y and V (x, y, a). That is, find Af (x,a), 
Af (y,a) and Af (V (x, y, a) ,a). 

Step 6. Find the rotation axis u = ( ui,u 2 ,u 3 ) where 


_ A (yiZ 2 ~ V 2 Zl) _ -A (xiz 2 + x 2 z 1 ) _ A (xiy 2 ~ X 2 Vl) 

Ul ai ■ Af (V (x, y, a) , 0 ) ’ a 2 • Af (V (x, y, a) , 0 ) ’ Ul a 3 • A/"(V (x,y, a) ,a) ' 

Step 7. Find the elliptical rotation angle using 

B(x,y,a ) 
cos t! = — — 

v / AA(x,a) v / AA(y,a) 

and define C = cos 0 and 5 = sin 0 where S' = \/l — C 2 . 

Step 8. Find the elliptical rotation matrix that rotates a; = (aq,yi,yi) to y = (x 2 , y 2 , z 2 ) 
elliptically on the ellipsoid aix 2 + a 2 y 2 + a 3 z 2 = 1 using Rodrigues matrix 


i? (u, C, S, A) 


aiu 2 + (l—aiu 2 ) C 

AusS <r i\ 

— a\u\u 2 (G-l) 


0 2 
Au 2 S 

03 


—CL1U1U3 (C-l) 


where u = (m, u 2 , u 3 ). 
Step 9. Define the matrix 


AU3S 

ai 


a 2 u\u 2 (C-l) 


a 2 tt|+ (l-a 2 u 2 ) C 


Am 5 

«3 


a 2 u 2 u 3 (C-l) 


A u 2 S 
a\ 

Au\S 

a 2 


a 3 uiu 3 (C-l) 
-a 3 u 2 u 3 (C-l) 


a 3 u 2 + (l -a 3 uf) C 


% (v,a) 


1 

AA (v,a) 


A/ - (v,a) — 2aiv 2 
— 201 ^ 2^1 
-2ain 3 ni 


-2a 2 uin 2 
A7 (v,a) — 2a 2 v\ 
-2 a 2 v 3 v 2 


—2a 3 v\v 3 
—2 a 3 v 2 v 3 
Af (v,a) — 2a 3 Ug 


for a given u = (vi,v 2 ,v 3 ). 

Step 10. Find the elliptical rotation matrix that rotates x = (xi,yi,yi) to y = (x 2 ,y 2 ,z 2 ) 
elliptically on the ellipsoid aix 2 + a 2 y 2 + a 3 z 2 = 1 using Householder matrices 


R (x, y,a)=U (y,a) H (x + y,a). 

Step 11. Define the set of Elliptic Quaternions H aija2!a3 = {q = r/o+g , ii+(/ 2 j+q , 3 k, qq. q\, q 2 , q 3 e 
M} with 


and 


i 2 = — a\ 

, j 2 = -02, 

k 2 = —a 3 

A 

A 

A 

—k = —ji, 

jk = — i = -kj, 

ki = — 

03 

01 

a 2 


Step 12. Find c = cos — = \j — r ^— and s = \/l — c 2 . Define the quaternion 


q = cos - + u sin - = c + s-iqi + su 2 j + su 3 k 

where 6 is the elliptical rotation angle and u = (ui,u 2 , u 3 ) = u\i + u 2 j + u 3 k is the rotation axis 
obtained in Step 6 and Step7. 

Step 13. Find the elliptical rotation matrix that rotates x to y elliptically on the ellipsoid 
using the matrix 
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R{q,a, A) 


qo+aiqj - a 2 ql ~ a 3 q f 
g 0 93 A 

2a\qiq 2 + 2- 

02 A 

o n 9092 A 

2ai(7i(?3 — 2- 


0 0 9o93 A 

2 a 2 q\q 2 — 2- 

ai 

9o - aiq\+a 2 ql - a 3 q% 

o . 0 9 o9iA 

2a 2 Q293 + 2- 


9092 A 

2a39i93 + 2- 

ai A 

o 0 9o9i A 

2a 3 g 2 93 - 2- 

a 2 

9o - °l9l - 0292+0393 


corresponding to q = qo + 9ii + 92j + 93 k. 

Remark 

In the n > 3 dimensional spaces, rotations can be classified such as simple, composite, and 
isoclinic, depending on plane of rotation. Simple rotation is a rotation with only one plane of 
rotation. In a simple rotation, there is a fixed plane. The rotation is said to take place about this 
plane. So points do not change their distance from this plane as they rotate. Orthogonal to this 
fixed plane is the plane of rotation. The rotation is said to take place in this plane. On the other 
hand, a rotation with two or more planes of rotation is called a composite rotation. The rotation 
can take place in each plane of rotation. These planes are orthogonal. In M 4 it is called a double 
rotation. A double rotation has two angles of rotation, one for each plane of rotation. The 
rotation is specified by giving the two planes and two non-zero angles (3 and 9 (if either angle is 
zero, then the rotation is simple). Finally, the isoclinic rotation is a special case of the composite 
rotation, when the two angles are equal [33] ■ In the 4 dimensional Euclidean and Lorentzian 
spaces, a skew symmetric matrix is decomposed as A = 9\A\ + 9 2 A 2 using two skew-symmetric 
matrices A\ and A 2 satisfying the properties A\A 2 = 0, A'\ = —A± and A% = — A 2 . Hence, the 
Rodrigues and Cayley rotation formulas can be used to generate 4 dimensional rotation matrices 

([33], m, HE], 0, and [J5]). 
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